Experimental evidence of a fractal dissipative regime in high-Tc superconductors 
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We report on our experimental evidence of a substantial geometrical ingredient characterizing 
the problem of incipient dissipation in high-Tc superconductors (HTS): high-resolution studies of 
differential resistance-current characteristics in absence of magnetic field enabled us to identify and 
quantify the fractal dissipative regime inside which the actual current-carrying medium is an object 
of fractal geometry. The discovery of a fractal regime proves the reality and consistency of critical- 
phenomena scenario as a model for dissipation in inhomogeneous and disordered HTS, gives the 
experimentally-based value of the relevant finite-size scaling exponent and offers some interesting 
new guidelines to the problem of pairing mechanisms in HTS. 
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Local inhomogeneities characterize HTS both on 
nanoscopic (e.g., periodic or aperiodic variation of 
local oxygen stoichiomctry) and mcsoscopic Q (e.g., oxy- 
gen depleted grain boundaries) spatial scales. While all 
the consequences on normal and superconducting charge 
transport in the former case has not been entirely clari- 
fied yet the case of grain boundaries is better understood: 
at least in a broad range of experimental parameters the 
supercurrent transport in polycrystalline samples relies 
on 'weak link network' (WLN), i.e., on mesoscopic su- 
perconducting islands interconnected by Josephson in- 
teraction. Although the transport features on nanoscopic 
scale may significantly differ from those characterizing a 
rather simple WLN problem (e.g., local inhomogeneities 
seem to give rise, as reviewed by Refs. [Q,^, to con- 
ducting stripes, clusters, wires or filaments which are, at 
least to some extent, mobile, compared to predominantly 
static weak links) the experimental evidence in favor of 
a-b plane Josephson junctions Q indicates that, besides 
qualitative similarities, the intrinsic and WLN transport 
are more closely related one to another than it had been 
foreseen earlier. 

Irrespective of the extent the processes at nano- and 
mesoscopic scales are related, the problem of charge 
transport in WLN represents an autonomous subject of 
much interest due to its relevance for general understand- 
ing of transport in heterogeneous systems and in Joseph- 
son Junction Arrays (JJA) in particular Focusing 
to the problem of dissipation there are convincing ar- 
guments, particularly in absence (or in small) magnetic 
fields, that the onset of dissipation is dominated rather by 
a phenomenon of percolation than the dynamical features 
of flux lattice [QJ^]. In a disordered-bonds (DB) model 
1^,^ the critical current Ic characterizing the dissipation 
onset reflects the connectivity threshold Pc of classical 
percolation networks ||9|,[Tot (such that Pc — pile)) so that 
the experimentally documented power-law-like current- 
voltage (I-V) characteristics can be naturally interpreted 
as a current-induced but in essence traditional critical 



phenomenon. Consequently, I-V characteristics should 
also reveal various manifestations of crossover between 
the relevant length scales known to underly the critical 
behavior. We show in this report that the latter crossover 
may be detected and quantitatively investigated in exper- 
imental I-V curves. In particular, we claim that the I-V 
characteristics are generically composed of the three dis- 
tinct regimes: a regime revealing no practical dissipation, 
a regime obeying conventional correlation length scaling 
(homogeneous regime) , and an intermediate regime obey- 
ing finite-size scaling (fractal regime). While the dissipa- 
tion in the former regimes has been already a subject to 
experimental reports and appropriate modeling |^,^ , the 
experimental results concerning the fractal regime are re- 
ported here for the first time. 

The basis of the model is the idea that the increasing 
current applied to disordered WLN decreases the fraction 
of Josephson-current-carrying bonds in a random man- 
ner. Hence, the applied current plays the role of random 
generator which in classical random electrical networks 
changes the relative fractions of their components. In DB 
model the elements and the relevance of this analogy 
has been studied in details. Here we focus to the problem 
of relevant length scales. In analogous classical networks 
there are two of them [^,0: the correlation length 
(the representative size of growing ramified clusters ) for 
p away from p^ and the sample size L for p close to p^ ( "at 
critical! ty" ) . In approaching p^, ^ diverges involving ex- 
ponent v (w 0.88 in 3D systems) and the power-law form 
of static and dynamic quantities manifests the property 
of spatial scaling. In particular, the resistance R of the 
random-superconductor network (RSN) disappears |ic| ] 
as (pc — p)^ oc , where s is the breakdown expo- 

nent (s « 0.8 in 3D). Close to pc the homogeneous-to- 
fractal transformation of the geometry of incipient clus- 
ter takes place and, while R becomes independent of p, 
the finite-size scaling relation U^M R{L) cx L~^/'^~^ re- 
places the ordinary R{L) cx one. In applying a sim- 
ilar scenario to current-induced transition in WLN we 
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FIG. 1. High-resolution difFerential resistance in Ag/HTS 
composite tape, core thickness « 30/im. The pronounced 
anomaly above the sharp dissipation onset at Ion is attributed 
to the sample size scaling > L). For higher currents the 
usual correlation length scaling (oc ^"^Z") takes over. Exper- 
imentally, dV/dl oc {I — Ic)^ , Ic = 0.6,2; = 2, as shown by the 
thin line. Vertical arrow illustrates the resolution: Its length 
corresponds to dc voltage of IbnV . The inset: V — I curve 
obtained by numerical integration of the measured dV/dl ~ I. 



assume that, in approaching 1^ (i.e., Pc) from above, the 
representative size of the largest phase coherent cluster 
diverges as ^ as well. The experimental studies [Q of 
the related homogeneous regime were shown to be in a 
close agreement with the predictions of the model. How- 
ever, the precise interpretation of the characteris- 
tic exponent (experimentally, dV/dl exponent is close to 
2) is still unresolved (see discussion section). Consid- 
ering the experimental accessibility of the crossover to 
the fractal regime we note that the unit of length in- 
volved in the WLN problem (i.e., average grain size £) 
belongs to mesoscopic (pm) scale. Hence, the observa- 
tion of a size effect introduced by competing length scales 
seems, for polycrystalline samples of reduced but still 
macroscopic thickness, as an open possibility. Indeed, 
the first high resolution dynamical resistance measure- 
ments |12j (achieving the equivalent voltage resolution of 
better than InV) on polycrystalline samples which are 
thin by their very design, HTS superconductor/normal 
metal composite tapes (superconducting core thickness 
in the range 10-50 fJ-m), revealed the two characteristic 
currents. As illustrated by Fig. 0, the lower one trig- 
gers the onset of low-level, non-exponential dissipation 
(onset current. Ion) while the higher one parameterizes 
the scaling behavior, dV/dl oc (/ — Ic)^,x « 2, of sub- 
sequently rapidly growing dissipation (therniodynamical 
critical current, /c). One could assume that a rather 
broad dissipative range between /o„ and Ic corresponds 
to validity of ^ > L when the incipient dissipative sites 



would fill the sample-sized network of fractal geometry. 
Indeed, the saturation-like behavior of dissipation in that 
range is, while in obvious disagreement with any flux- 
creep model (exponential in applied current), at least 
in qualitative agreement with general independence of 
any observable [dV/dl in our case) on ^ in the range of 
sample-sized fractal [^jlO|. A similar observation of the 
broad range of low-level dissipation in composite tapes 
has also been reported by other authors but interpreted 
by less fundamental causes p^ . 

By studying a thickness dependence in appropriate 
samples we prove now the presence of geometrical con- 
straints of fractal nature in initial dissipation of HTS 
in a more quantitative way (the results on composite 
tapes. Fig. represent just a qualitative indication). 
The presence of many spurious and/or overlapping ef- 
fects in composite tapes and HTS films precludes 
obtaining a firm quantitative information on critical be- 
havior from these samples. We performed therefore 
the measurements of I-dV/dl characteristics on a well- 
characterized, non-textured (i.e., isotropic) polycrys- 
talhne RBaaCuaOr-^r (R=Y,Gd) bulk sample (a WLN 
prototype!) in many successive steps, after its thickness 
had been gradually reduced by fine plan-parallel grind- 
ing. In that way, apart from various thickness, all the 
measurements were performed on the same initial sam- 
ple. The transport properties of the sample (e.g., room 
temperature resistivity, Tc, resistive transition width) 
did not change in all stages of its thickness. The mea- 
surements we report on in this paper covered the sample 
thickness range of 20-1000 iim (factor of 50). For thick- 
nesses above approximately 60 /im only the unique 'ther- 
modynamic' critical current, accompanied by the usual 
power-law-like (oc (/ — Ic)^ ,x k, 2) growth of dissipation 
(specific for scaling regime in a very large sample, ^ <C -L), 
have been detected. In the sample stages involving all 
smaller thicknesses the two characteristic currents. Ion 
and Ic, have been observed, just as in composite tapes. 
Some of the experimental dV/dI{L) curves were shown in 
Fig. 1^ using moderate (main figure) and a very high dy- 
namical resistance resolution (inset) . The anomalous dis- 
sipative range between Ion and Ic is rather complex but 
systematically depends on sample thickness: The size of 
the onset anomaly drastically increases by decreasing L. 
The analysis and interpretation have been performed in- 
side the DB model |^,^ which provides both the limiting 
behavior in the correlation length scaling range ^ L 
(thin lines in Figs. |], ||) and the estimate of the width of 
the range of sample size scaling (^ > L). The crossover 
between these two ranges takes place when the diverging 

C(p) — ^ Pc\^^ , becomes equal or higher than the 
sample size. The unit of length is ^, the network unit cell 
size (for i we used i ~ 5fim, the average grain size of the 
sample). In other words, as long as the fraction of 'good' 
bonds deviates from pc (percolation threshold of an infi- 
nite system) by Ap = {i/Ly/" = (£/L)i i36 or less, the 
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FIG. 2. Thickness dependence of differential resistances 
of polycrystalline GdBa2Cu307_i: sample (open symbols) 
and homogeneous-to- fractal phase boundary (thick grey line). 
The thickness (in ^m) is designated by numbers near the 
experimental curves. The thin lines represent the predicted 
power- laws characterizing the homogeneous regime (see text). 
A sizable deviation which scale the sample size corresponds 
to sample-sized cluster of fractal geometry. Large triangles 
were used to extract the Rc{L) scaling in Inset to Fig. 3. In- 
set: Temperature dependence of the very onset of anomalous 
dissipation in the thinnest (19.7/im) sample. 



macroscopic properties should have a weak dependence 
on p (i.e., on current in our case) and the underlying 
ramified sample-sized cluster should be, geometrically, a 
fractal. 

In DB model j4|J|] a linear p{I) approximation has 
been shovirn to work well close to Pc (but still outside 
Ap): Pc- p ^ {c2/ci){pc/Ic){I ~ Ic) where C2/C1 (= g 
henceforth) is a geometrical factor of order 1. The cur- 
rent interval compatible with Ap, A/, reads therefore 
A/ = Icg{Ap/pc) = (/c5/pc)(^/i)^-"^- The onset cur- 
rent can be now defined as Ion — Ic — A/ and the corre- 
sponding current density is expected to depend strongly 
on sample thickness L : Jon — >/c(l — {g/Pc){^/ L)^'^^^)- 
Experimentally, while the determination of /o„ in thin 
samples is quite straightforward (the onset of dynamical 
resistance is very sharp, inset to Fig. ||) the determination 
of Ic is not; Ic represent just a parameter in DB model for 
dissipation Q], dV/dl = Rf{gpc/IcT{I - IcY , where Rf 
represents the total resistance of WLN in homogeneous 
regime. 

An interesting observation is that the thermodynamic 
critical current density Jc is strongly thickness depen- 
dent as well. Fig. |. The presence of a size-effect in 
Jc is, however, a rather well-known [p^ although some- 
what neglected phenomenon. We found that the ob- 
served increase of Jc for a factor of 5 by thinning the 
cross-section fits perfectly the general breakdown formula 
Jc oc 1/(1 + {alogioLY) which predicts that the 
critical current density vanishes in thermodynamic limit. 
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FIG. 3. Thickness dependence of thermodynamical critical 
current (Jc) and onset current (Jon) density- experimental 
points (symbols) and model predications (thin lines). Inset: 
Thickness dependence of dV/ dl at Ic (at criticality) indicates 
the fractal geometry of the network. The slope for homoge- 
neous (Euclidean) electrical networks is n = 1. 



Our results are compatible with Jc{L) — 5Jc{Lmax)/{^ + 
{logio{L/Lmin)Y) where L„iax and L.nm are maximal 
and minimal sample thickness, respectively. Inserting the 
latter expression into the derived one for Jon (compris- 
ing only one adjusting parameter, g) one gets the peaked 
curve. Fig. ^ as a prediction of this model. The two 
dependencies (for Jc and Jon) joins smoothly in increas- 
ing L as Ap (i.e.. A/) continuously vanishes, as well as ^ 
itself, in this limit. The remarkable overlap of experimen- 
tal points and the model predictions, Figs. 2, 3, illustrate 
the reality of the model, in spite of its simplicity. It is 
interesting to note that Fig. ^ can be interpreted as a 
kind of phase diagram: the thick grey line separates the 
homogeneous from the fractal phase of the cluster inside 
which the initial dissipation grows. 

The main quantitative results of this work is plotted 
as an Inset to Fig. ||. It plots the value of dynamical 
resistance in fractal phase Rc{Ic) for each available sam- 
ple thickness L. The particular current at which Rc has 
been taken was Ic, the representative of Pc in current- 
induced transitions. The Rc{L) relationship can be fit 
nicely by a power- law Rc{L) oc L~". The exponent value 
(n — 3.36) deviates strongly from the value n = I char- 
acterizing homogeneous networks (the scaling with n = 1 
is strictly obeyed in, e.g., Rf{L) dependence). Also, the 
quantity vin — 1) which in finite-size scaling calculations 
gives the dynamical exponent of homogeneous regime is 
numerically very close (= 2.1) to experimentally well- 
documented Q value X — 2 valid in that regime. Both 
arguments provide therefore the evidence for fractal ge- 
ometry involved in initial dissipation. 
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The important issue which remains to be clarified is the 
interpretation of the value of the exponent n, as well as of 
the related exponent x characterizing the homogeneous 
regime of I-V characteristics. The experimental values 
of exponents (a; ~ 2, n w 3.4) are, while mutually consis- 
tent, in clear disagreement with those obtained by identi- 
fying RSN and WLN, i.e., 0.8, n = s/v+l « 1.9. 
There could be several reasons why the WLN exponent, 
X, may differ from the classical one, s. A well-known 
example is the 'Swiss-cheese' morphology in continuum 
percolation [Q, equivalent to the case of broad distri- 
bution in bond resistances, shown to influence the ex- 
ponent. The other is experimentally documented non- 
universal conduction in carbon-black-polymer composite 
1^ attributed to peculiarities of tunneling as a mech- 
anism of local conduction. Both the broad distribution 
and tunneling seems as natural possible causes for expo- 
nent deviation in WLN of HTS. We also note that there 
are some obvious differences between current-generated 
(WLN) and random-generated (RSN) clusters. Better 
understanding of these differences could probably come 
from very recent and exciting studies of self-organized 
pof and 'small- world' |^ networks. 

The observation of a fractal dissipative regime offers 
an interesting new guideline towards understanding the 
intrinsic pairing interaction in HTS. There are namely 
numerous arguments that the intrinsic intra- and inter- 
plane charge transport takes place actually in a hetero- 
geneous conductive medium 0|, with percolation play- 
ing probably the important role as well |Q. Moreover, 
these conditions are considered, according to some au- 
thors |2^, as substantial ingredients of the mechanism 
of superconductivity itself. The involved heterogeneity 
may rely either on charge separation, stripes, wires, etc., 
I ^ cluster formation 1^ or on filamentary fragmentation 
1 22[ . Under these circumstances it seems quite reasonable 
to assume that the intrinsic current transfer may include 
the fractal network as well (at least in certain range of 
relevant transport parameters). Given the electrically 
heterogeneous local properties, combined with vicinity 
of metal-insulator transition, the associated elastical (vi- 
brational) network (which is formally isomorphic to its 
electrical counterpart |l^) might be not only heteroge- 
neous but may posses a fractal geometry as well. In 
a fractal elastical lattice the vibrations are, instead of 
extended phonons, the localized high-frequency fractons 
p3| , p4| which may contribute in pairing. Relaying on pe- 
culiarities of the fracton density of state ||2^ , such as high 
cut-off frequency and/or high-frequency 'missing modes' 
psf , the pairing temperatures could be higher than those 
associated to classical phonons. On basis of our results 
we suggest therefore consideration of a fractal dynami- 
cal lattice as a possible source of non-standard pairing 
interactions at high temperatures. 
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